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ABSTRACT 



Context. The chemical evolution of galaxies is investigated within the framework of the star formation rate (SFR) dependent integrated 
galactic initial mass function (IGIMF). 

Aims. We study how the global chemical evolution of a galaxy and in particular how [a/Fe] abundance ratios are affected by the 
predicted steepening of the IGIMF with decreasing SFR. 

Methods. We use analytical and semi-analytical calculations to evaluate the mass-weighted and luminosity-weighted [a/Fe] ratios in 
early-type galaxies of different masses. 

Results. The models with the variable IGIMF produce a [a/Fe] vs. velocity dispersion relation which has the same slope as the 
observations of massive galaxies, irrespective of the model parameters, provided that the star formation duration inversely correlates 
with the mass of the galaxy (downsizing). These models also produce steeper [a/Fe] vs. cr relations in low-mass early-type galaxies 
and this trend is consistent with the observations. Constant IMF models are able to reproduce the [a/Fe] ratios in large elliptical 
galaxies as well, but they do not predict this change of slope for small galaxies. In order to obtain the best fit between our results 
and the observations, the downsizing effect (i.e. the shorter duration of the star formation in larger galaxies) must be milder than 
previously thought. 

Key words. Stars: abundances - stars: luminosity function, mass function - supernovae: general - Galaxies: evolution - Galaxies: 
elliptical and lenticular, cD - Galaxies: star clusters 



1. Introduction 

It is nowadays widely accepted that most stars in galaxies form 
in star clusters (Tutukov I19781 Lada & Lada 2003). This has 
been observed in a number of different galaxies; from the Milky 
Way to the dwarf galaxies of the Local Group (Wyse et al. |2002| 
Massey I2003I Piskunov et al. I20041 ). Within each star cluster, 
the initial mass function (IMF) can be well approximated by the 
canonical two-part power-law form £(m) oc mT a (e.g. Pflamm- 
Altenburg Weidn er & Kroupa l2007l hereafter PWK07). Massey 
& Hunter (I1998I I have shown that for stellar masses m > a few 
M a slope similar to the Salpeter ( 119551 index (i.e. a = 2.35) 
can approximate well the IMF in clusters and OB associations 
for a wide range of metallicities, whereas many studies have 
shown that the IMF flattens out below m ~ 1 M (Kroupa, Tout 
& Gilmore lT993l ChabrierEOQD. 

On the other hand, star clusters are also apparently dis- 
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tributed according to a single-slope power law, £ eC ] 
where M ec i is the stellar mass of the embedded star cluster. There 
is a general consensus that this slope /3 should be of the order 
of ~ 2 (Zhang & FaU [19991 Lada & Lada |2003t Hunter et al. 
2003 ), although a /3 as high as 2.4 can also be realistic (Weidner, 
Kroupa & Larsen 2004). According to this correlation, small em- 
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bedded clusters are more numerous in galaxies. They provide 
therefore most of the stars but not most of the massive ones, 
since they are preferentially formed in massive clusters (Weidner 
& Kroupa 2006]). As a consequence of this mass distribution of 
embedded clusters, the integrated IMF in galaxies, the IGIMF, 
can be steeper than the stellar IMF within each single star clus- 
ter (Kroupa & Weidner[2003] Weidner & Kroupa|^05l). 

The Salpeter IMF slope has been used in a very wide range 
of modelling, providing good fits with observations concern- 
ing the cosmic star formation history (Calura, Matteucci & 
Menci 2004), the X-ray properties of elliptical galaxies (Pipino 
et al. 120051 ). the chemical evolution of dwarf galaxies (Larsen, 
Sommer-Larsen & Pagel 2001) and of the Milky Way (Pilyugin 
& Edmunds 119961 but see also Romano et al. 2005). Broadly 
speaking, a flatter than Salpeter IMF produces a larger frac- 
tion of massive stars. The large production of oxygen (and 
of a-elements in general) leads to lower [Z/O] metallicity ra- 
tios. A steep IMF slope would instead be biased towards low- 
and intermediate-mass stars, underproducing oxygen and there- 
fore resulting in larger [N/O] and [C/O] abundance ratios. On 
the other hand, iron will also be overproduced compared to a- 
elements, since it comes mainly from Type la SNe which orig- 
inate from C-O deflagration of binary systems of intermediate 
mass. Therefore, galaxies characterized by a steep IMF will tend 
to have [a/Fe] ratios lower than models in which the IMF is flat. 

The scenario of a variable integrated galactic initial mass 
function (IGIMF) has been applied in models of chemical evolu- 
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tion (Koppen, Weidner & Kroupa 2007 '), producing an excellent 
agreement with the mass-metallicity relation found by Tremonti 
et al. d20041 > . However, these authors consider only the effect of 
the IGIMF on the global metallicity and the evolution of abun- 
dance ratios has not yet been explored in the literature. In a se- 
ries of papers we plan to study the impact of the IGIMF on the 
abundance ratios in different classes of galaxies, using different 
methodologies. In this paper we study, by means of simple ana- 
lytical and semi-analytical models, the evolution of [a/Fe] ratios 
in galaxies, in particular in early-type ones. It is, in fact, now well 
established that the [a/Fe] ratios in the cores of elliptical galaxies 
increase with galactic mass (Weiss, Peletier & Matteucci 1995; 
Kuntschner et al. 1200 11 1 and this poses serious problems to the 
current paradigm of hierarchical build-up of galaxies (see e.g. 
Thomas et al. l20031 here after T HOM05; Nagashima et al. l2005l 
Pipino, Silk & Matteucci 2009 1 Calura & Menci, in preparation). 
In fact, in the classical hierarchical models the most massive el- 
lipticals take a longer time to assemble and therefore form stars 
for a longer time than less massive galaxies, thus producing a a 
trend of [a/Fe] vs. mass which is opposite of what is observed 
(see Thomas, Maraston & Bender 2002; Matteucci 120071 1 . 

We will show in this paper that the trend of increasing [a/Fe] 
vs. galaxy mass is naturally accounted for in models of elliptical 
galaxies in which the IGIMF is implemented. The second paper 
of this series will be devoted to the study of the chemical evolu- 
tion of the Solar Neighborhood and of the local dwarf galaxies 
and in this case we will make use of detailed chemical evolution 
models. Another paper of this series will study the evolution of 
galaxies by means of chemodynamical models, in order to ana- 
lyze how the IGIMF changes the feedback of the ongoing star 
formation in galaxies and how this affects the chemical evolu- 
tion. 

The plan of the present paper is as follows. In Sect. 2 we 
summarize the IGIMF theory and the formulations we adopt. In 
Sect. 3 we describe how we calculate the Type la and Type II SN 
rates in galaxies in which the SFR is given. Once we know the 
Type la and Type II SN rates, it is possible to calculate the [a/Fe] 
ratios. This has been done in Sect. 4 for ellipticals and early-type 
galaxies in general. A discussion and the main conclusions are 
presented in Sect. 5. 

2. The determination of the integrated galactic 
initial mass function 

The determination of the IGIMF has been described previously 
(Kroupa & Weidner |2003l Weidner & Kroupa |2005l PWK07). 
The IGIMF theory is based on the assumption that all the stars 
in a galaxy form in star clusters. Surveys of star-formation in 
the local Milky Way disk have shown that 70 to 90 % of all 
stars appear to form in embedded clusters (Lada & Lada 2003 1 
Evans et al. 2008 ). The remaining 10-30 % of the apparently dis- 
tributed population may stem from a large number of short-lived 
small clusters that evolve rapidly by dissolving through energy 
equipartition and residual gas expulsion. It is therefore reason- 
able to assume that star formation occurs in embedded clusters 
with masses ranging from a few M upwards. The IGIMF, inte- 
grated over the whole population of embedded clusters forming 
in a galaxy, becomes 

/~M cc ,, m „(SFR(r)) 

&HMF(»t; SFR(f)) = f(m < mmax^eclCMecOdAfeci, (1) 

J McCain 

where M ec i. m i n and M ec i, max (SFR(f)) are the minimum and max- 
imum possible masses of the embedded clusters in a population 



of clusters and m max = m max (M ec i) (eqs.|3]and|4]i. For M ec i, min we 
take 5 M (the mass of a Taurus-Auriga aggregate, which is ar- 
guably the smallest star-forming "cluster" known), whereas the 
upper mass of the embedded cluster population depends on the 
SFR. The correlation between M ec i >max and SFR has been deter- 
mined observationally (Larsen & Richtler 2000; Weidner et al. 
2004) and can be expressed with the correlation 

log Med^im = log^ML + 0.75 log^r + 6.77, (2) 

where tfr is the SFR in M yr~' and &ml is the mass-to-light 
ratio, typically 0.0114 for young stellar populations (Smith & 
Gallagher 2001). This empirical finding can be understood to 
result from the sampling of clusters from the embedded cluster 
mass function given the amount of gas mass being turned into 
stars per unit time (Weidner et al. 2004). 

The stellar IMF (i.e. the IMF within each embedded cluster) 
has the canonical form %(m) = km~ a , with a = 1 .3 for 0.08 M < 
m < 0.5 M and a = 2.35 (i.e. the Salpeter slope) for 0.5 M < 
m < m max , where m max depends on the mass of the embedded 
cluster. In order to determine m max and the proportionality con- 
stant k we have to solve the following two equations (Kroupa & 
Weidner |2003l : 

/-("max 

M ec \ = I m^(m)dm, (3) 

Jmi ov 

I f(m)dm = l, (4) 

where m\ ovl is the smallest considered stellar mass (0.08 M in 
our case) and m max » is the upper physical stellar mass and its 
value is assumed to be 150 M (Weidner & Kroupa l2004l l. Eq.|4] 
indicates that, by definition of m max , there is only one and exactly 
one star in the embedded cluster with mass M ec \ whose mass is 
larger than or equal to m max . 

The last ingredient we need is the distribution function of 
embedded clusters, £eci(Meci)> which, as we have mentioned in 
the Introduction, we can assume proportional to . In this 
work we have assumed 3 possible values of j3: 1.00 (model 
BETA 100), 2.00 (model BETA200) and 2.35 (model BETA235). 
In Fig. [T]we have plotted the resulting IGIMFs for different val- 
ues of SFR. In particular, we have tested 20 SFRs, ranging from 
10~ 4 to 100 M yr _1 , equally spaced in logarithm. To appreciate 
better the differences between various models, we have plotted in 
Fig.|2]IGIMFs for 3 different values of the SFR: SFR 10" 2 M 
yr" 1 (heavy lines), SFR 1 M yr" 1 (middle lines), SFR * 10 2 
M yr (light lines). We have considered all the possible values 
of p. model BETA 100 (dashed lines), BETA200 (dotted lines), 
BETA235 (solid lines). For clarity, we have plotted the IGIMFs 
only for masses larger than ~ 2 M , since in the range of low 
mass stars the IGIMFs do not vary. As expected, the model with 
the steepest distribution of embedded clusters (model BETA235) 
produces also the steepest IGIMFs. This is due to the fact that 
model BETA235 is biased towards embedded clusters of low 
mass, therefore the probability of finding high mass stars in this 
cluster population is lower. 

We can also notice from Fig. |2]that the differences between 
IGIMFs with SFR 1 M yr 1 (middle lines) and SFR * 10 2 
M yr~' (light lines) are not very pronounced. This is due to the 
fact that for both these SFRs, the maximum possible mass of the 
embedded cluster is very high (see eq.|2]i, therefore in both cases 
the upper possible stellar mass of the whole galaxy is very close 
to the theoretical limit of 150 M . This can be seen in Fig. [3] 
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Fig. 1. IGIMFs for different distributions of embedded clusters 
(e.g. different values of ft). Upper panel: model BETA235; cen- 
tral panel: model BETA200; lower panel: model BETA 100. For 
each panel we have considered 20 possible values of SFR, rang- 
ing from 10~ 4 M yr -1 (lowermost curves) to 100 M yr _I (up- 
permost curves), equally spaced in logarithm. 



(lower panel) in which we plot the variation of ra max as a func- 
tion of SFR as deduced from eqs. [4] and [3] This correlation is 
valid for all the possible values of ft because it is determined by 
£(m) and not by % ec i(M ec i). As we can see from Figs. Q] and [2] the 
IGIMFs are characterized by a nearly uniform decline, which 
follows approximately a power law, and a sharp cutoff when m 
gets close to m max . In Fig. [3](upper panel) we plot therefore also 
the slope that better approximates the IGIMF in the range 3-16 
M . This is the range of masses where most of the progenitors of 
SNell and SNela originate (see Sect. 3). Of course, the steeper 
the distribution of embedded clusters is, the steeper the corre- 
sponding IGIMFs are. Fig. [3] shows also what we have noticed 
before, namely that the various IGIMFs saturate for SFR > 1 M 
yr . Finally, in Fig. [3] (middle panel) k a is shown as a function 
of the SFR for the various models. k a is the number of stars per 
unit mass in one stellar generation (see e.g. Greggio 2005) and 
its value is given by 



kfy 



m ^ 1G1MF ( m )dm 



(5) 



This parameter is useful to calculate the SNII rates (see Sect. 3). 

3. The determination of Type la and Type II SN rates 

3.1. Type II SN rates 

Stars in the range m up < m < m max (where m up is the mass limit 
for the formation of a degenerate C-O core) are generally sup- 
posed to end their lives as core-collapse SNe. These SNe divide 
into SNell, SNelb and SNelc according to their spectra. For our 
purposes, this distinction is not useful and we will suppose that 



Fig. 2. IGIMFs for different SFRs: 
' M yr 

We have considered all the possible values of ft: model BETA 100 
(dashed lines), BETA200 (dotted lines), BETA235 (solid lines). 



10~ 2 M yr _1 (heavy 
lines); 1 M yr -1 (middle lines); 10 2 M yr -1 (light lines). 




-2 
log (SFR/M s yr->) 



Fig. 3. Lower panel: m max (in M ) as a function of SFR (in M 
yr -1 ). Middle panel: k a (number of stars per unit mass in one 
stellar generation). Upper panel: IGIMF slopes in the range 3 - 
16 M . Notations as in Fig. [2] 



all the core-collapse supernovae are indeed SNell. These SNe 
produce the bulk of a-elements and some iron (one third approx- 
imately). The standard value of m ap is 8 M but stellar mod- 
els with overshooting predict lower values (e.g. Marigo 1200 H . 
However, stars more massive than m up can still develop a de- 
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generate O-Ne core and end their lives as electron-capture SNe 
(Siess 2007). We will assume for simplicity that all the stars with 
masses larger than 8 M end their lives as SNell, therefore the 
SNII rate is simply given by the rate at which massive stars die, 
namely: 



Rs 



tit) 



-f 



flit - r m )^GiMF[m, if/it - T m )]dm, 



(6) 



where if/ is the SFR and t,„ is the lifetime of a star of mass m. 
The lifetime function is adapted from the work of Padovani & 
Matteucci ( fT993l , 



f 1.2m- 1 - 85 + 0.003 Gyr if m > 6.6 M G 
\l0 /(m) Gyr ifm<6.6M , 



where 



Sim) 



[0.334- VI -79 - 0.2232 x (7.764 - log(m)) 



0.1116 



(7) 



(8) 



In eq. [6] the IGIMF is calculated by considering the SFR at the 
time t - t,„, therefore it depends both on time and on mass. 

It is instructive to analyze models in which SFR is constant 
during the whole evolution of the galaxy. In this way, eq.|6] sim- 
plifies into 



R 



SNII 



f 



£igimf('«, ft)dm = if/k a n >& 



(9) 



where, as we have seen in Sect. 2, k a is the number of stars per 
unit mass in one stellar generation (therefore tf/k a is the number 
of stars formed per unit time) and n>8 is the number fraction 
of stars with masses larger than 8 M . In Fig. [4] (lower panel) 
we plot Rsnii (in cen 1 ) as a function of SFR for the 3 adopted 
values of ft. It is worth pointing out that, even if SFR is constant, 
Rs nii starts increasing only after the star with mass m max ends 
its life and reaches a constant value only after the lifetime of 
a 8 M star. This lifetime (~ 28 Myr with the adopted lifetime 
function) is however negligible compared to the Hubble time, 
therefore it is reasonable to consider Rsnii constant with time. 
For SFRs larger than ~ 10~ 2 M yr _1 Rsnii increases almost 
monotonically with SFR, whereas it drops dramatically for SFR 
< 10~ 2 M yr~'. This is mostly due to the drop of n>g at low 
SFRs (Fig. |H upper panel) which in turn depends on the fact 
that the upper mass m max for these values of SFR starts reducing 
significantly and it gets very close to 8 M for a SFR of 10~ 4 
M yr~' (Fig.O, therefore only a very narrow interval of stellar 
masses gives rise to SNII explosions. From Fig. [3] we can see 
instead that the variation of k a with SFR is not very significant, 
therefore k a affects only mildly the Type II SN rates. 

It is nowadays getting popular to consider SN rates normal- 
ized to the stellar mass of the considered galaxy. The usually 
chosen unit of measure is the SNuM (1 SNuM = 1 SN cen -1 
10 -10 M» -1 , where M* is the current stellar mass of the galaxy). 
In this case, models in which the SFR is constant cannot attain 
a constant Type II SN rate in SNuM since the stellar mass of 
the galaxy increases with time. We therefore calculated Rsnii in 
SNuM as a function of time for the various models. The stellar 
mass of the galaxy at each time t is given by tfi- f <m ( t) dt, where 
f<m(t) is the mass fraction of stars, born until the time t that have 
not yet died. 

Fig. [5] shows the evolution with time of the Type II SN rate 
for different models and different SFRs, assuming a constant 
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Fig. 4. Lower panel: Rsnii (in cen~') as a function of SFR (in 
M yr _1 ). Upper panel: «>8 (number fraction of stars more mas- 
sive than 8 M ) as a function of SFR. Notation as in Fig. [2] 



SFR for 14 Gyr. These results are compared with the average 
SNell rates (in SNuM), observationally derived by Mannucci et 
al. (|2()05) in SOa/b galaxies (solid boxes), Sbc/d galaxies (dotted 
boxes) and irregular ones (dashed boxes). We can notice that, 
for the models BETA 100 and BETA200 only the mildest SFRs 
can reproduce the final SNII rates in SOa/b galaxies, whereas 
model BETA235 can fit the final SNII rate of SOa/b galaxies for 
a wide range of SFRs. On the other hand, all the models predict 
final SNII rates significantly below the observations of irregu- 
lar galaxies and the final values for model BETA235 fail also to 
fit the observed rates in Sbc/d galaxies. It is important to note, 
however, that the stellar mass in galaxies is usually calculated as- 
suming some (constant) IMF. Under the assumption that the IMF 
changes with the SFR, the determinations of the stellar masses 
must be revisited. PWK07 showed that the IGIMF effect (i.e. the 
suppression of the number of massive stars with respect of low- 
mass stars) can be very significant in dwarf galaxies, whereas 
in large galaxies it tends to be very small. Moreover, a constant 
SFR for 14 Gyr is not a reasonable description of the star forma- 
tion history of irregular (and Sbc/d) galaxies which often expe- 
rience an increase of the SFR in the last Gyrs of their evolution 
(see e.g. Calura & Matteucci 2006). For this reason, the calcu- 
lated SNII rates of late type galaxies tend to fit the observations 
at smaller ages. 

3.2. Type la SN rates 

In order to calculate the SNIa rates, we assume the so-called 
Single Degenerate Scenario of SNIa formation. It is commonly 
assumed that a SNIa explodes when a C-O white dwarf in a bi- 
nary system reaches the Chandrasekhar mass after mass accre- 
tion from a companion star. According to the Single Degenerate 
channel of SNIa explosion, the accretion of matter occurs via 
mass transfer from a non-degenerate companion (a red giant or a 
main sequence star) filling its Roche lobe (Whelan & Iben |1973l . 
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to reach the Chandrasekhar mass. With these assumptions, the 
limits of integration in eq. ( fTOl t are: 



5 10 
t (Gyr) 



5 10 
t (Gyr) 



5 10 
t (Gyr) 



Fig. 5. Rsnii (in SNuM) as a function of time (in Gyr) for mod- 
els with constant SFR and different values of /?: model BETA 100 
(left panel), BETA200 (central panel), BETA235 (rigth panel). 
In each panel, the lowermost curve has SFR = 10~ 4 M yr _1 , 
the uppermost has SFR = 10 2 M yr~' and the other SFRs are 
equally spaced in logarithm. Observations of the core-collapse 
SN rates in SOa/b galaxies (solid boxes), in Sbc/d galaxies (dot- 
ted boxes) and in irregular ones (short-dashed boxes) are also 
shown. Data are taken from Mannucci et al. (20Q5J 

In this way, the SNIa rate depends on the number distribution of 
C-O white dwarfs, but also on the mass ratio between primary 
and secondary stars in a binary system. The SNIa rate in the 
framework of the Single Degenerate Scenario has been analyt- 
ically calculated by a number of authors assuming a universal 
IMF (see Valiante 2009 and references therein). Here we follow 
the formulation of Greggio & Renzini (1983) and Matteucci & 
Recchi d2001| l but we modify it to take into account that, in the 
framework of the IGIMF, the IMF changes according to the SFR. 
The SNIa rate in this case turns out to be: 



fWB.inf = max[2m2(t), 3M ] 
m B ,su P = 8M + m 2 (t), 



// m in = max 



RiM 



J /( 



f(lJ.)^(t-T m2 )^i G]MF [m B , if/(t-T„ l2 )]d[i dm B , (10) 



where A is a normalization constant (assumed to be 0.09 in the 
following). Although theoretical arguments demonstrate that A 
should be small (e.g. Maoz 2008) its value is usually calibrated 
with the Milky Way. Unfortunately, our analytical approach does 
not allow us to simulate the Milky Way within the IGIMF theory, 
therefore we take 0.09 as a reference value and postpone a more 
careful discussion about it to the follow-up numerical paper (but 
see also Sec. [4] for a study of the variation of A for early-type 
galaxies). m B is the total mass of the binary system, ni2 is the 
mass of the secondary star, yu = mil m B and f{fi) is the distribu- 
tion function of mass ratios (see below). It is commonly assumed 
that the maximum stellar mass able to produce a degenerate C- 
O white dwarf is a 8 M star, therefore the maximum possible 
binary mass is 16 M . The minimum possible binary mass is as- 
sumed to be 3 M in order to ensure that the smallest possible 
white dwarf can accrete enough mass from the secondary star 



m 2 (t) m B - 8M G 



m B 



(11) 
(12) 



(13) 



The distribution function of mass ratios is generally described 
as a power law (/(//) // y ), but the value of y is still much 
debated in the literature (Duquennoy & Mayor fl~99 1 1 Shatsky 
& To kovinin 120021 Kouwenhoven et al. 120051 and therefore we 
will take it as a free parameter. 

Fig. [6] shows the evolution with time of the Type la SN rate 
for different models and different SFRs, analogously to Fig.[5]for 
SNell rates. Also shown (dashed lines) for comparison are SNIa 
rates obtained for a model with fixed (i.e. not SFR-dependent) 
IMF. We assume the canonical stellar IMF (i.e. the IMF within 
each embedded cluster) which, as mentioned in Sect. [2] has the 
form £(m) = km'", with a = 1.3 for 0.08 M < m < 0.5 M Q 
and a = 2.35 above 0.5 M . As we can see, at large SFRs 
model BETA 100 produces rates almost indistinguishable from 
the ones obtained with the fixed canonical IMF (see also Kroupa 
& Weidner l2003l l. In this figure y is assumed to be 2 (Tutukov 
& Yungelson 1980). This large value of y favors the occurrence 
of SNela in binary systems with similar masses. Such a steep 
mass ratio distribution that favors equal-mass binaries may result 
from dynamical evolution of stellar populations in long-lived 
star clusters (Shara & Hurley |2002| l. We can notice again that 
only model BETA235 at very low SFRs seems able to reproduce 
the SNIa rates in SOa/b and E/S0 galaxies. However, we point 
out that the comparison with the observed SNIa rates in ellip- 
tical galaxies is meaningless because they have stopped form- 
ing stars several Gyr ago, therefore they have evolved passively 
since then. For them we cannot therefore assume a constant SFR 
for 14 Gyr (see Sect. EJ. On the other hand, model BETA235 
produces SNIa rates that only match the observed rates in dwarf 
irregular galaxies at their peak. Therefore, assuming y — 2, the 
best value for /3 seems to be 2 (but see the comment in Sect. 13.11 
about the possible inconsistency of the published determination 
of stellar masses, at least for irregular galaxies). To show the de- 
pendence of the results on y we show in Fig. |7]the SNela rates 
obtained assuming y = 0.3. This flatter distribution function im- 
plies that a larger fraction of binary systems with small mass 
ratios end up as SNela. We can notice from this figure that the 
observed SNIa rates in spiral galaxies are reproduced by a larger 
range of SFRs, whereas the disagreement with the observed rates 
in irregular galaxies worsens. 



4. A test of the IGIMF: [a/Fe] ratios in early-type 
galaxies 

The study of the average stellar [a/Fe] ratio in galaxies repre- 
sents an important constraint for our models, since this quan- 
tity depends both on the adopted galactic star formation history 
and on the stellar IMF (Matteucci 2001). In local ellipticals, the 
observed correlation between the central velocity dispersion <x, 
which reflects the total stellar mass, and the stellar [a/Fe] is inter- 
preted as due to the shorter star formation timescales in the most 



1 Note that usually observative papers adopt mass ratios q = mi/nii 
instead of fi (Abt & Levy 1976). The mass ratio distribution function 
f(q) can be obtained from f(p.) by means of a simple change of variable. 
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Fig. 6. RsNia (in SNuM) as a function of time (in Gyr) for mod- 
els with constant SFR. Notations and symbols as in Fig. [5] with 
the addition of the observed SNIa rate in E/SO galaxies (long- 
dashed boxes) and the SNIa rates predicted assuming a constant 
IMF (dashed lines). Data are taken from Mannucci et al. (2005). 
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Fig. 7. As in Fig.|6]but for y = 0.3. 



massive galaxies (Pipino & Mattuecci 2004; THOM05) which in 
turn implies also that the most massive galaxies experience the 
most intense episodes of star formation. For this reason, the aver- 
age stellar [ar/Fe] vs <x relation represents a valuable test for the 
IGIMF, since the IGIMF is a function of the galactic star forma- 
tion rate. The issue of a variable IMF among elliptical galaxies to 
explain the [a/Fe] vs. <x relation has already been explored with 



Fig. 8. Duration At of the star formation (assumed constant) as a 
function of SFR (lower scale) and cr (upper scale) assuming the 
Af-luminous mass relation of THOM05 (see text). 



success by Matteucci (1994) but assuming ad hoc variations of 
the IMF slope. In this section we test, using well-established and 
observationally constrained star formation histories of early-type 
galaxies of various masses, if the physically motivated IGIMF 
can equally well reproduce this correlation. 

To simplify the calculations, the SFR is assumed to be con- 
stant over a period of time At. We have numerically tested that 
this crude approximation about the star formation history does 
not affect significantly the results. The value of At as a function 
of galaxy luminous mass is adopted from the work of THOM05, 
who, on the basis of the observational relation between [a/Fe] 
and cr, showed the existence of a downsizing pattern for ellipti- 
cal galaxies, according to which the smaller ellipticals form over 
longer timescales (see also Matteucci 1994; Cowie et al. 1996; 
Kodama et al. 2004). Since the present-day stellar mass is given 

rAt 

in this case by ij/ ■ J fi ow (t)dt (where /i ow (0 is the fraction of 
long-living stars, namely the stars, born at the time t, that live 
until the present day), it is possible to derive a relation between 
the SFR and the duration of the star formation activity At, which 
we show in Fig. [8] This relation saturates at 14 Gyr since this is 
assumed to be the age of the Universe. A similar relation can be 
recovered from the work of Pipino & Matteucci (2004) assuming 
that the star formation occurs only until the onset of the galac- 
tic wind, however the two SFR-Af relations do not significantly 
differ. 

For each galaxy (characterized by a specific SFR over a pe- 
riod At) we calculate the average yield from SNell of a chemical 
element i, 



y? 



J^" lm " x yi(m)^i G i MF (m, f)dm 



(14) 



where yi(m) is the yield of chemical element i produced by a sin- 
gle star of mass m. The nucleosynthetic prescriptions are taken 
from Woosley & Weaver dl995l l. We have however halved the 
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Fig. 9. IGIMF- averaged SNII yields of oxygen (solid lines), 
iron (dotted lines) and magnesium (dashed lines) as a function 
of SFR (in M yr _1 ) for different values of /?: model BETA 100 
(left panel), BETA200 (central panel), BETA235 (right panel). 

iron yields, in accordance with Timmes, Woosley & Weaver 
( 1995) and Chiappini, Matteucci & Gratton (119971 1. because it is 
known that only in this way it is possible to reproduce the [a/Fe] 
in Galactic stars. Unfortunately, Woosley & Weaver (1995) cal- 
culated yields only for stellar masses up to 40 M . We assume 
that the yields of stars more massive than 40 M are equal to 
the 40 M star yields. Given the very limited amount of stars 
in the range 40 M Q < m < m max the results are not sensitive to 
this assumption. In Fig. [9] we show the IGIMF-averaged SNII 
yields of oxygen (solid lines), iron (dotted lines) and magne- 
sium (dashed lines) as a function of the SFR for different values 
of (3. For what concerns SNela, we assume the yields reported by 
Gibson, Loewenstein & Mushotzky d!997| l, based on the work of 
Thielemann, Nomoto & Hashimoto (1993). 

Once we know the SNIa yields and the IGIMF-averaged 
SNII yields for each galaxy, we can calculate the mass fraction 
jr(f) (where a is O or Mg) produced until the time t by using 
the formula: 



(15) 



where Ri d and /?snii are the SNIa and SNII rates, given by eqs. 
[I0]and|9] respectively, and y Ia are the SNIa yields. 

At this point, we can compute the theoretical average stellar 
abundances by means of: 



[a/Fe] = logj 



J ■'At 
o jr e {i) ■ Aow(i)dt 



M„ 



l0gl °W 



(16) 



where M to , is the total present-day stellar mass and a and Fe Q 
are the solar abundances of a-elements and Fe, respectively, 
taken from Anders & Grevesse ( 1989). This means that our the- 
oretical abundance ratios represent values mass-averaged over 



the stars which survive until the present time (see also Thomas, 
Greggio & Bender [1999j ). 

The observable in elliptical galaxies is the velocity disper- 
sion instead of the mass, so in order to properly compare our 
results with observations we need to assume a correlation be- 
tween the stellar mass and the velocity dispersion of galaxies 
(Faber-Jackson relation). We assume: 



0.86M 



0.22 



(17) 



(Burstein et al. 1997), where cr is the velocity dispersion in km/s. 
The resulting relation between cr and At can be seen from Fig. [8] 
where we have indicated in the upper scale the cr corresponding 
to each SFR. 

In Fig. [10] we show our results for y = 0.3 comparing our 
models with observations taken from THOM05 and references 
therein (filled squares). We can first notice that, as expected, 
the model BETA 100 (heavy dashed lines), giving rise to flat- 
ter IGIMFs (see Fig. [JJ(, produces larger [a/Fe] ratios. In fact, 
flatter IGIMFs result in a larger fraction of massive stars and, 
therefore, to a larger production of a-elements. We can also ap- 
preciate that the models reproduce quite well the [a/Fe] (both 
[O/Fe] and [Mg/Fe]) ratios in elliptical galaxies, at least for the 
models BETA200 and BETA235. To appreciate the effect of the 
IGIMF approach, we plot also (long-dashed line) a model with 
the fixed canonical IMF which, as mentioned in Sect. [2] and re- 
minded in Sect. 13.21 has the form g(m) = krrT a , with a = 1.3 
for 0.08 M < m < 0.5 M and a = 2.35 above 0.5 M . The 
curves obtained with the IGIMF tend to flatten out at large cr, 
whereas the curve obtained with the constant IMF shows a con- 
stant slope. This demonstrates once more that the adoption of the 
IGIMF is particularly remarkable in the low-mass (and low-c) 
galaxies. The curve with a constant IMF approaches asymptot- 
ically the model BETA 100 since this model at large SFRs pro- 
duces the flattest IMFs (see Fig. [JJ). Besides a small shift of a 
few tenths of dex (which can be fixed increasing the parame- 
ter A in Eq. QJJi, the curve with a constant IMF reproduces well 
the trend of [a/Fe] vs. cr of the THOM05 sample, demonstrating 
that the downsizing (or inverse- wind) models (Matteucci 1994; 
Pipino & Matteucci 2004) are also perfectly capable of explain- 
ing this trend in large elliptical galaxies. However, evidence is 
mounting that [a/Fe] ratios in early-type dwarf galaxies are so- 
lar or sub-solar. For instance, van Zee, Barton & Skillman (2004) 
showed that [a/Fe] ratios (derived from Lick indices) of a sample 
of Virgo dwarf irregular galaxies range between -0.3 and solar. 
Also in the cluster Abell 496 the smallest galaxies show [Mg/Fe] 
to be solar or sub-solar (Chilingarian et al. 2008). To show that 
we have also plotted in Fig. [10] (open triangles) the data of a 
sample of low-mass early-type galaxies by Sansom & Northeast 
(2008). These data confirm that the [a/Fe] vs. cr relation is prob- 
ably steeper in the low-mass regime and that our IGIMF results 
can naturally predict this behavior. However, in order to properly 
test our results in the low-mass regime more data are needed. 

It is worth pointing out that in this figure (and in the follow- 
ing ones) we have considered only model galaxies for which the 
SFR is smaller than 100 M yr~' . This is the reason why the data 
points reach larger cr than the results of our model. In extreme 
starbursts the IMF might become top-heavy as evident by the 
mass-to-light ratios in ultra-compact dwarf galaxies which are 
ultra-massive "star clusters" that form when the SFR is very high 
(Dabringhausen, Kroupa & Baumgardt 2009 ) and this will need 
to be incorporated in the IGIMF calculations (work in prepara- 
tion). 

In general, in local early-type galaxies the stellar abundances 
are measured by means of various absorption-line Lick indices, 
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log cr(km/s) 



Fig. 10. Mass- weighted [Mg/Fe] (upper panel) and [O/Fe] (lower panel) vs. cr for models with a constant SFR over a period of time 
At which depends on the stellar mass of the galaxy (see text) assuming y = 0.3. Heavy solid line: model BETA235; heavy dotted 
line: model BETA200; heavy short-dashed line: model BETA100; long-dashed line: model with a constant (i.e. not SFR-dependent) 
canonical IMF. The filled squares are the observational values and relative error-bars as reported by THOM05 and references therein 
and the light solid line is the least-square fit of this data. The open triangles are [ff/Fe] ratios reported by Sansom & Northeast ( 2008 ) 
and the light short-dashed line is the corresponding least-square fit. 



such as Mg b and < Fe >= 0.5(Fe52720 + Fe5335) (THOM05). 
To properly compare predictions to observational abundance 
data obtained for local ellipticals, in general one should derive 
the luminosity-weighted average abundances. The real abun- 
dances averaged by mass are larger than the luminosity-averaged 
ones, owing to the fact that, at constant age, metal-poor stars 
are brighter (Greggio 1 1997b . To calculate the luminosities we 
have made use of the Starburst99 package (Leitherer et al. 1999, 
Vazquez & Leitherer 2005), producing L(t) for each value of 
SFR and /3. The results are shown in Fig. [TT|for the first 100 
Myr (the luminosities remain almost constant after 100 Myr). 



As expected, since model BETA 100 is characterized by the flat- 
test IGIMFs, it produces also the largest luminosities. We have 
then calculated luminosity-weighted mass ratios by using the 
formula: 

t f- e (t)L(t)dt aQ 

[a/Fe] = log 10 l ogl0 (18) 

JE Ut)dt Fe ° 

The results are shown in Fig.[T2] As we can see, the results dif- 
fer very little (by a few hundredths of dex at most) compared 
to the mass-averaged abundance ratios. We have checked these 
results also using the spectro-photometric code of Jimenez et 
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Fig. 11. Stellar luminosities (in L Q ) as a function of time 
(in Myr) for models with constant SFR and different values 
of p. model BETA 100 (left panel), BETA200 (central panel), 
BETA235 (right panel). In each panel, the lowermost curve has 
SFR = 10~ 4 M yr" 1 , the uppermost has SFR = 10 2 M yr" 1 and 
the other SFRs are equally spaced in logarithm. These luminosi- 
ties have been obtained with the Starburst99 package (Leitherer 
et al. [19991 Vazquez & Leitherer|2005]l 



al. (119981 1 (see also Calura & Matteucci 120031) but the results 
do not differ appreciably compared with the ones obtained with 
the Starburst99 package. Indeed, it has been already shown in 
the literature (but for constant IMFs) that the discrepancy for 
the [Mg/Fe] ratio computed by averaging by mass and by lumi- 
nosity is very small, with typical values of 0.01 dex (Matteucci, 
Ponzone & Gibson 1998; Thomas et al. 1999). We have con- 
firmed this finding also in the case of the IGIMF. 

To check how much our results depend on the assumption of 
a variable Af with stellar mass, we plot in Fig.[13]the [a/Fe] ob- 
tained assuming a constant value of Af = 1 Gyr. The agreement 
with the observations is still quite good; in particular the mod- 
els maintain an increasing trend of [ar/Fe] with cr. However, the 
curves tend to flatten out too much at larger cr, at variance with 
the trend shown by the observations. This is due to the fact that, 
as pointed out in Sect. [2] the various IGIMFs ferrates of star for- 
mation larger than 1 M yr~' do not show very large differences. 
Therefore, the assumption of a star formation duration inversely 
proportional to the stellar mass of the galaxy (or in other words 
the downsizing) is a key ingredient to understand the chemical 
properties of large elliptical galaxies. 

To appreciate the dependence on the distribution function of 
mass ratios in binary stars (the parameter y introduced in Sect. 
|3.2| > we plot in Figs. [14] and Q3] the results of models with y = 
2.0 and y = -0.3, respectively. The curves obtained with y = 2.0 
tend to be slightly steeper than the ones shown in Fig. [10] (and 
slightly steeper than the observations) but the agreement remains 
still good, in particular for the models BETA 100 and BETA200. 
If we assume y = -0.3 an excellent match with the observations 
is instead provided by the model BETA235. Models BETA 100 



Fig. 12. As in Fig. QJjJbut with luminosity-weighted [a/Fe] ra- 
tios, calculated by means of eq. [18] 




1.4 1.6 1.8 2 2.2 2.4 

log cr(km/s) 



Fig. 13. As in Fig.[l0]but with Af = 1 Gyr for each stellar mass. 



and BETA200 show the same slope of the observational data but 
shifted by a few tenths of dex. A slight increase of the parameter 
A in eq. [10] would make these models perfectly compatible with 
the observations. 

It is particularly remarkable that the trend of [a/Fe] ratios 
vs. cr (namely an increase of [a/Fe] with cr) is naturally repro- 
duced using the IGIMF approach, without any further assump- 
tion or fine-tuning of parameters. This is for instance at variance 
with what hierarchical clustering models of structure formation 
would tend to produce, since in this case larger elliptical galaxies 
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Fig. 14. As in Fig.QUbut with y = 2.0. 
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Fig. 15. As in Fig.[T0]but with y = -0.3. 



are formed later, out of building blocks in which the [a/Fe] ratio 
has already dropped (e.g. Thomas et aL 12002) . It is worth men- 
tioning that De Lucia et al. (2006 ), by means of a semi-analytical 
model adopting the concordance A CDM cosmology, suggested 
that more massive ellipticals should have shorter star formation 
timescales, but lower assembly (by dry mergers) redshift than 
less luminous systems. This is one of the first works based on the 
hierarchical paradigm for galaxy formation producing downsiz- 
ing in the star formation histories of early-type galaxies through 
the inclusion of AGN feedback (see also Bower et al. 2006; 
Cattaneo et al. 2006), although they did not compute the [afFe]- 



Fig. 16. Normalization constant A to adopt in order to obtain 
the best fit with the observational data as a function of y for 
model BETA 100 (heavy dashed line), BETA200 (heavy dotted 
line) and BETA235 (heavy solid line). The light lines represent, 
for each value of f3, the normalized chi square of the represented 
model (scale on the right axis). 



cr relation for ellipticals. However, the lower assembly redshift 
for the most massive system is still in contrast to what is con- 
cluded by Cimatti, Daddi & Renzini (2006), who show that the 
downsizing trend should be extended also to the mass assem- 
bly, in the sense that the most massive ellipticals should have 
assembled before the less massive ones. Very recently, Pipino 
et al. (120081 ) showed that even in semi-analytical models able 
to account for the downsizing, the [a/Fe] vs. cr relation is not 
reproduced. 

Although the agreement between our results and the obser- 
vations is good, none of the models presented so far fits per- 
fectly the data at low and high cr simultaneously. In order to 
work out an overall best model, for each value of y and f3 we 
have checked, by means of a minimization of the normalized chi 
square, which normalization constant A fits better the data. The 
results are shown in Fig. [16] As we can see, model BETA235 
seems to be preferable and the best agreement between data and 
models is obtained for the model BETA235 with y = 2.5 and A 
= 0.036. In general, the best fits are obtained with large values of 
y, although that requires low values of A. A large value of y, fa- 
voring equal-mass binary systems, is consistent with the results 
of Shara & Hurley (2002), although observational surveys cited 
in Sect. [32] seem to indicate lower values of y. 

We should however not forget that the Af-luminous mass re- 
lation we have used in this work has been obtained by THOM05 
assuming a constant IMF. We have therefore checked, starting 
from our best model, namely a model with fi = 2.35, y = 2.5 
and A = 0.036, how this relation should change in order to best 
fit the data. It turns out that, within the IGIMF theory, the best 
Af-luminous mass relation is given by: 



logAf = 2.38 -0.24 log M„_ 



(19) 
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Fig. 17. Mass-weighted [Mg/Fe] (upper panel) and [O/Fe] 
(lower panel) vs. cr for our best IGIMF model (namely model 
BETA235 with y = 2.5 and A = 0.036) with a Af-luminous 
mass relation described by eq.[l9](heavy solid lines). Also shown 
(heavy long-dashed lines) a model with fixed, canonical IMF, 
y — 2.5 and A = 0.13. Notations and symbols as in Fig.fTUl 

where At is in Myr and M tot in M Q , and this relation produces 
the results shown in Fig. [17] We have also plotted in this figure 
a model with canonical stellar IMF (heavy long-dashed line) in 
which we have increased the value of A up to 0.13 in order to 
better reproduce the data. We can notice again that the canonical 
IMF can perfectly fit the [a/Fe] ratios in large ellipticals but it 
shows a constant slope and therefore it cannot equally well re- 
produce the [a/Fe] ratios in dwarf galaxies, at variance with the 
IGIMF model. The comparison between relation[T9land eq. 5 of 
THOM05 is displayed in Fig. [18] We can notice here that the 
downsizing effect (namely the shorter duration of the star for- 
mation in larger galaxies) is milder, in the sense that the At for 
large galaxies is (slightly) larger than the timescale calculated 
by THOM05, whereas at low cr the star formation durations are 
significantly lower than the ones predicted by THOM05. 

5. Discussion and conclusions 

In this paper we have studied, by means of analytical and semi- 
analytical calculations, the evolution of [a/Fe] ratios in early- 
type galaxies and in particular their dependence on the lumi- 
nous mass (or equivalently on the velocity dispersion cr). We 
have applied the so-called integrated galactic initial mass func- 
tion (IGIMF; Kroupa & Weidner |2U03l Weidner & Kroupa 
2005) theory, namely we have assumed that the IMF depends 
on the star formation rate (SFR) of the galaxy, in the sense 
that the larger the SFR is, the flatter is the resulting slope of 
the IGIMF. This kind of behavior would naturally tend to form 
more massive stars (and therefore more SNell) in large galax- 
ies, which are characterized by more intense star formation 
episodes. Therefore, it is expected that, since a-elements are 
mostly formed by SNell, the most massive galaxies are also the 



Fig. 18. Af-luminous mass relation obtained with eq. [19] (solid 
line) and derived by THOM05 (their eq. 5; dashed line). 



ones which attain the largest [a/Fe] ratios, in agreement with the 
observations. One of the main aims of this paper was to quanti- 
tatively check whether the chemical evolution of galaxies within 
the IGIMF theory is able to accurately fit the observed [a/Fe] vs. 
cr relation. 

We have analytically calculated the SNII and SNIa rates with 
the IGIMF assuming 3 possible slopes of the distribution func- 
tion of embedded clusters, £ ec ] oc Afjj, where M ec i is the stellar 
mass of the embedded star cluster; in particular we have consid- 
ered £ = 1.00 (model BETA 100); = 2.00 (model BETA200); 
/? = 2.35 (model BETA235). We have seen that, if we consider 
constant SFRs over the whole Hubble time, the final SNIa and 
SNII rates agree quite well with the observations of spiral galax- 
ies (in particular the SOa/b ones). The agreement with the ob- 
served rates in irregular galaxies is not good, but a constant SFR 
over the whole Hubble time is not likely in irregular galaxies, 
which probably have experienced an increase of the SFR in the 
last Gyrs of their evolution. 

To calculate the [a/Fe] ratios with the IGIMF we assumed 
that early-type galaxies form stars at a constant rate over a pe- 
riod of time At which depends on the total luminous mass of 
the considered galaxy. This hypothesis is based on the work of 
THOM05 who, on the basis of observational grounds, showed 
the existence of a downsizing pattern for elliptical galaxies, i.e. 
that the most massive galaxies are the ones with the shortest Af. 
We have then calculated the production of a-elements and Fe by 
SNell (in particular we have calculated IGIMF-averaged SNII 
yields) and by SNela and we have calculated mass-weighted and 
luminosity-weighted [a/Fe] ratios for each model galaxy, char- 
acterized by different SFRs and [}. 

The resulting mass-averaged [a/Fe] vs. cr relations show the 
same slope as the observations in massive galaxies as reported 
by THOM05, irrespective of the value of (3 and of the distribu- 
tion function of mass ratios in binaries f(fi) oc ^ (which affects 
the SNIa rates), although models with [3 = 2.35 and large val- 
ues of y seem to be preferable. Some models show a shift (of 
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a few tenths of dex) compared with the observations but this 
can be fixed increasing (or decreasing) the fraction, A, of binary 
systems giving rise to SNela, which is an almost unconstrained 
parameter. It is however remarkable that all the models we have 
calculated show the same trend of the observations because if, as 
commonly argued, large elliptical galaxies form out of mergers 
of smaller sub-structures (hierarchical clustering), it would be 
natural to expect that they are the ones with the lowest [a/Fe] ra- 
tios because they form later, out of building blocks where [a/Fe] 
has already dropped. 

It is worth pointing out that the [cr/Fe] ratios do not depend 
on the gas flows (infall and outflow) experienced by the galaxy 
(Recchi et al. 2008) therefore our results do not depend on spe- 
cific infall and outflow parameters, which make them particu- 
larly robust. However, these parameters affect the overall metal- 
licity of the galaxy, therefore they need to be taken into account 
in order to check whether our models can correctly reproduce 
the mass-metallicity relation. As mentioned in the Introduction, 
Koppen et al. (2007) have already shown that the IGIMF the- 
ory is able to reproduce the mass-metallicity relation found by 
Tremonti et al. (2004) in star-forming galaxies. We are checking, 
by means of detailed numerical models, that the IGIMF theory is 
able to reproduce at the same time the mass-metallicity relation 
and the [c/Fe] -cr relation in early-type galaxies. This study will 
be presented in a forthcoming paper. 

We have also considered models in which the IMF does not 
vary with the SFR and, because of the variations of Af with SFR, 
these models are compatible with the observations of large ellip- 
tical galaxies as well. However, these models produce a [a/Fe] 
vs. cr relation that can be described as a single-slope power- 
law, whereas the IGIMF models bend down significantly at low 
masses (and low cr). This is because the IGIMF becomes partic- 
ularly steep in the galaxies with the mildest SFRs and this adds 
to the downsizing effect (namely the decreasing duration of the 
SFR with increasing mass). From our study therefore, an impor- 
tant conclusion is that a very reliable observable to test the valid- 
ity of the IGIMF theory is the observation of the [ar/Fe] ratios in 
dwarf galaxies. The available data on [o7Fe] ratios in low-mass 
early-type galaxies show indeed some steepening of the [o7Fe] 
vs. cr relation, in agreement with the IGIMF predictions. 

We have also tested how much our results depend on the as- 
sumption of a variable Af with stellar mass by computing mod- 
els with Af = 1 Gyr irrespective of the stellar mass. The agree- 
ment between models and observational data is still reasonably 
good but the curves tend to flatten out too much at large stel- 
lar masses compared with the observations (and with the IGIMF 
models). This indicates that the downsizing remains a fundamen- 
tal ingredient to understand the chemical properties of early-type 
galaxies. However, if we check for which Af-luminous mass re- 
lation we obtain the best fit between data and models, it turns 
out that the downsizing effect must be milder than predicted 
by THOM05, in the sense that large galaxies form stars for a 
slightly longer timescale than calculated by THOM05, whereas 
low-mass galaxies have star formation durations significantly 
shorter. Although the exact form of the best-fit Af-luminous 
mass relation is subject to a number of parameters (IGIMF pa- 
rameters; parameters regulating the SNIa rate etc.) and might 
change once larger and more detailed abundance measurements 
are available, the result of a milder downsizing effect compared 
to the findings of THOM05 is robust. 

Eventually, we have seen that luminosity-weighted [a/Fe] ra- 
tios agree very well with the mass-weighted ones (with relative 
differences of a few hundredths of dex at most), in accordance 
with the results of Matteucci et al. (I1998I I. 



We remind the reader that, with our analytical approach to 
chemical evolution, we are making some important simplifying 
assumptions. For instance, our computation of the interstellar 
jr(f) given by eq. [151 does not take into account in detail the 
lifetimes of massive stars. Furthermore, our present calculations 
do not take into account the variation with time of the metallicity 
in galaxies, which should also influence the stellar yields. From 
the various tests performed so far, and from the comparison of 
our results with numerical results (Thomas et al. 1999; Pipino & 
Matteucci 2004), we have verified that these assumptions may 
play only some minor role in determining the zero-point, but not 
the slope of the predicted [a/Fe] vs cr relation. All of these sim- 
plifying assumptions will be relaxed in our forthcoming paper, 
where we will present a numerical approach to the role of the 
IGIMF in galactic chemical evolution. 

The main results of our paper can be summarized as follows: 

- Models in which the IGIMF theory is implemented naturally 
reproduce an increasing trend of [a/Fe] with luminous mass 
(or cr), as observed in early-type galaxies. 

- However, models with constant duration of the star forma- 
tion produce a [a/Fe] vs. cr relation which flattens out too 
much at large cr. Only models in which the star formation 
duration inversely correlates with the galactic luminous mass 
(downsizing) can quantitatively reproduce the observations. 

- Models in which the IGIMF is implemented show (at vari- 
ance with the constant IMF models) a steepening of the 
[a/Fe] vs. cr relation for small galaxies, therefore the IGIMF 
theory can be tested by observing the [a/Fe] in dwarf galax- 
ies. The observations available so far are in agreement with 
our predictions. 

- Luminosity-weighted abundance ratios differ from the mass- 
weighted ones by a few hundredths of dex at most. This re- 
sult, already known for constant IMF models, has been con- 
firmed in the IGIMF framework. 

- In order to obtain the best fit between our results and the 
observed [a/Fe] ratios in early-type galaxies, the downsizing 
effect (namely the shorter duration of the star formation in 
larger galaxies) has to be milder than previously thought. 

- The best results are obtained for a cluster mass function 
£ ed oc M~^ 35 , indicating that the embedded cluster mass 
function should have a Salpeter slope. 
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